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for r - 

h(d)=d*(1/(200*pi))*asin(sqrt(2)*d/(4000)); 
z(d)=2*d* 6.99/(2828) -h(d) ; 
hi(d)=d*(1/(200*pi))*asin(d/(4000)) ; 
z2(d)=2*d*10/(4000)-hi(d); 

h2(d)=d*40/4000-h(d); 
zl(d)=2*d*30/(4000)-h2(d); 
h3(d)=d* 60/4000-h(d); 
z3(d)=2*d*50/4000-h3(d); 
h4(d)=d*80/4000-h(d); 
z4(d)=2*d*70/4000-h4(d); 
end 

for d=l:4500 

yl(d)=10*d/4000; 
y2(d)=30*d/4000; 
y3(d)=50*d/4000; 
y4(d)=70*d/4000; 
end 

x=1:4500; 
d=l:4000; 

plot (d,h(d) ,d, z (d) ,d,hl (d) , d, z2 (d) ,d,h2 (d) , d, zl (d) ,d,h3 (d) ,d, z3 (d 
) ,d,h4 (d) ,d,h5(d) , ) 
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